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The problem of detecting an enemy submarine whose possible position was 
revealed by the hit of a torpedo is known as the “Flaming Datum” problem. All previous 
studies devoted to this theme make unrealistic assumptions about the speed of the 
escaping target when dealing with a diesel-electric submarine. In this kind of submarine 
the constraint imposed by the remaining charge of its batteries determines that its 
behavior is essentially conservative in how fast it should escape. 

The objective of this thesis is to explore the idea of varying speed in the flaming 
datum problem. Two different approaches are considered. An analytical model is 
developed based on the relationship among some of the physical factors that could 
determine or constrain the behavior of a diesel submarine while escaping from the area of 
the flaming datum. The second approach considers a discrete event simulation using the 
Java-based Simkit package. Data analysis is used to determine a possible fit for the 
simulation results. Several tactics are explored to determine their effects on detection 


probability. 
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EXECUTUVE SUMMARY 


The problem of detecting an enemy submarine whose possible position has been 
revealed by the hit of a torpedo on a surface vessel 15 known as the “Flaming Datum" 
problem. Since World War II, several studies have been devoted to the issue where a 
submarine tries to escape from the datum area while being hunted by a helicopter 
equipped with variable depth sonar (VDS) and/or antisubmarine weapons. All these 
studies make assumptions about the speed of the escaping target that are unrealistic when 
dealing with a diesel-electric submarine. The constraint imposed by the remaining 
charge of its battery determines that the behavior of this kind of submarine is essentially 
conservative in how fast it should escape. Even though the submarine is able to maintain 
a constant speed for long periods of time, this is only valid for very low speeds. Low 
speeds are not tactically acceptable for a submarine trying to evade a searcher, especially 
in the first minutes after the attack when operating in the vicinity of the Datum. 

According to opinions and experiences collected by the author, a diesel-electric 
submarine will tend to escape, as a general rule, as soon as possible from the zone where 
an attack recently occurred. This tactic is motivated by the efforts of the enemy to find the 
submarine, but it is only sustainable for short periods of time depending on the 
characteristics of the battery, and its remaining charge. 

The objective of this thesis is to develop and test a model for the Flaming Datum 
problem with a submarine trying to evade at varying speeds. First, the author develops 
an analytical model based on the existing relationship among some of the physical factors 
that determine or constrain the speed and behavior of a diesel submarine trying to escape 


from the area of the flaming datum. This model predicts the probability of detecting the 


Xi 


submarine as a function of the helicopter's delays in arriving at the datum area. To 
obtain a more accurate result for the same problem a discrete event simulation using the 
Java based Simkit package is developed and validated as believable. This simulation is 
also used to test some variations where either the submarine or the helicopter tries a 
specific tactic. 

Even though, both the analytical model and the simulation. show some 
weaknesses that may affect the credibility of the results, the possible sources of 
inaccuracy are identified and explained. It turns out that the submarine speed should be a 
strong function of time, rather than constant as usually supposed. 

This study represents a more realistic method of posing the Flaming Datum 
problem when dealing with diesel-electric submarines that must respect an energy 
constraint. This constraint was never considered in any other study. For this reason the 


developed models may be considered a valid starting point for future works. 
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І. INTRODUCTION 


A DESCRIPTION OF THE PROBLEM 


The problem of detecting an enemy submarine whose possible position has been 
revealed by the hit of a torpedo on a surface vessel is known as the “Flaming Datum” 
problem. Since World War II, several studies have been devoted to the issue where a 
submarine tries to escape from the datum area while being hunted by a helicopter 
equipped with variable depth sonar (VDS) and/or antisubmarine weapons. All these 
studies make assumptions about the speed of the escaping target that are unrealistic when 
dealing with a diesel-electric submarine. The constraint imposed by the remaining charge 
of its battery determines that the behavior of this kind of submarine is essentially 
conservative in how fast it should escape. Even though the submarine is able to maintain a 
constant speed for long periods of time, this is only valid for very low speeds. Low speeds 
are not tactically acceptable for a submarine trying to evade a searcher, especially in the 
first minutes after the attack, when operating in the vicinity of the Datum. 

According to opinions and experiences collected by the author, a diesel-electric 
submarine will tend to escape, as a general rule, as soon as possible from the zone where 
an attack recently occurred. This tactic is motivated by the efforts of the enemy to find the 
submarine, but it is only sustainable for short periods of time depending on the 
characteristics of the battery, and its remaining charge. 

The objective of this thesis is to develop and test a model for the Flaming Datum 

Problem with a submarine trying to evade at varying speeds. To do so, the author 


develops an analytical model based on the existing relationship among some of the 


physical factors that determine or constrain the speed and behavior of a diesel submarine 
trying to escape from the area of the flaming datum. The model is tested with a discrete 


event simulation. 


B. BACKGROUND 


The previous studies can be divided into those that have been used to solve the 
problem as either one-dimensional or two-dimensional games. In the first group, 


Meinardi [Ref. 1] considers a multi-staged search game where the submarine hides in a 
row Of Mo boxes. At each stage, the helicopter chooses a box and inspects it. If the target 


is not found, it can move on to an adjacent box or remain in the same one. Both target 
and searcher keep track of the searched boxes. The time spent to transit among the boxes 
is not taken into consideration, and thus, the applications of this game are limited for 
further analysis. 

Baston and Bostock [Ref. 2] develop a game that could be considered a 
continuous version of the flaming datum problem in one-dimensional space. The 
submarine moves in a long narrow channel and the helicopter searches for it by carrying 
a limited number of bombs. Each bomb has a lethal radius R and a limited number of 
bombs are available. These two characteristics can be considered analogous to the sonar 
detection range and the number of dips available. In this game the submarine and the 
helicopter are restricted in movement by their maximum speeds (U and V) and 
instantaneous changes in velocity are allowed. 

The second group deals with two-dimensional cases. Danskin [Ref. 3] supposes 


that the datum is known precisely and the submarine strategy is to choose a course 0 and 


a speed v with 0«- v«- vo, where vo is the maximum speed. It holds to both throughout 
the search and this establishes the submarine's position within its speed circle [Figure 1]. 
Danskin also supposes that the submarine does not know what the helicopter is doing, so 
that it does not need to alter its speed or course. The helicopter dips are assumed to be 
inside the submarine's speed circle and non-overlapping. Danskin's assumptions seem to 
be very drastic. 

Cheong [Ref. 4] develops a game where the submarine and helicopter are 
assumed to be moving at the constant speeds of U and V respectively. However, the 
submarine is allowed to execute evasive maneuvers by moving directly away from the 
last helicopter dip as well as moving away from the datum. Meanwhile the helicopter 


conducts a random search. 


Figure 1. Danskin approach. Submarine speed circle, showing a choice 
of strategy. (From Ref. 3, p. 511) 


Washburn [Ref. 5, p 2-7] states that expected number of detections for a random 


search in an expanding disk can be modeled as: 





| 2RV di (Equation 1) 
ay(t) 


г 


where И 15 Ше searcher's speed, R is the range of its cookie cutter sensor, t is the delay 
of the helicopter arriving at the datum, and y(t) is the expanding disk radius at time г, 


which is dependent on the submarine’s maximum speed. Equation 1 will be used in 
developing an analytical model in Chapter II. 

Washburn and Thomas [Ref. 6] model the flaming datum problem as a dynamic 
search game. They relax the rules of motion for the submarine, so that its future position 
could be anywhere inside its speed circle and not necessanly connected to the previous 
position. This assumption, in a sense, permits an unbounded target speed. In fact, the 
maximum speed is only considered to generate the radius of the area to be searched for, 
or the maximum speed circle. Even though this assumption appears to favor the 
submarine, Washburn and Thomas are of the opinion that it is of minor consequence. The 
authors compare the similarity of their results with Cheong’s, and state that this similarity 
might imply that both models are actually wrong. According to Washburn “no Flaming 
Datum Problem has ever been solved or even usefully bounded, so the required standard 


of comparison is missing” [Ref 6. P 421]. 


С. АРРКОАСН 


Our first approach is to use an analytical model to develop a plausible speed 
profile for the submarine. This model is based on the existing relationship among some of 
the physical factors that determine or constrain the speed and behavior of a diesel 
submarine trying to escape from the area of the flaming datum. Parameters and variables 
such as power, energy, speed, speed sustainability and time are employed to determine an 
expression to compute the probability of detection of the submarine as a function of time. 
The submarine's foe 1s assumed to be a generic helicopter equipped with variable depth 
sonar (VDS) with a “cookie-cutter” sensor. This helicopter arrives at the flaming datum 
after some delay 7 and has a limited endurance. If the submarine is not detected after 
this period, it will be considered to have escaped. 

In a second and more accurate approach a Monte Carlo discrete event simulation 
using the Java-based Simkit Software package is developed. A submarine with a set of 
continuous decreasing velocity functions that are consistent with the capabilities of a 
generic submarine used in the analytical model is considered. 

A detailed description of the assumptions and strategies for the submarine and 


helicopter are given in the next chapters. 
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Il. ANALYTICAL MODEL 


A. ASSUMPTIONS 


As mentioned in Chapter I, all previous studies make assumptions about the 
speed of the target that are unrealistic when dealing with a diesel-electric submarine. The 
constraint imposed on this kind of submarine by the remaining charge of the batteries 
forces its behavior to be essentially conservative when considering escape speed. We 


make the following assumptions: 


1. Submarine Assumptions 


a. The submarine will normally try to escape as quickly as possible in the 
first few minutes ın order to make the FOC (farthest on circle) larger, and 
to increase the probability of escape. This behavior is justified since the 
first reaction of an attacked force will be to verify that the submarine is no 
longer in the vicinity of the datum. 

b. The battery is initially close to having a maximum charge. This 1s a 
reasonable assumption for a deliberate attack. 

c. The submarine is able to estimate the expected helicopter’s delay in 
arriving in the datum area. 


d. The depth factor will be ignored. 


Helicopter Assumptions 


The position of the datum 1s known with no error. 

The searcher knows the submarine s capacities related to its speed. This 
assumption is considered valid when taking into account that open sources 
and commercial manuals are readily available. As a general rule, the 
searcher should know who the possible enemies are. 

The helicopter moves at a constant speed when not dipping. The times to 
dip and raise the sonar, and search in a position, are also constant. 

The helicopter does not assume any particular behavior for the submarine; 
therefore its dips are uniformly distributed in the submarine’s maximum 


possible distance circle (MPD). 


B. THEORETICAL FOUNDATION 


Rydill and Burcher [Ref 7, p 291], state that the effective power (PE) needed to 
"tow" the submarine at speed u(t) in a submerged condition 15 given by: 
1 
where Sper 15 the wetted area of the bare hull and casing plus the wetted area of all 


appendages minus the area covered by the appendages, C7, is the total resistance 


coefficient or drag coefficient and p is the sea water density. 


Now let 
1 
J = 5 PS rer Crs 


Thus, the energy of the submarine’s batteries required to move a submarine with a 


speed u(t) during an interval [0, £] can be expressed as: 


(у= (шау (Equation 3) 


Then, the constraint imposed by the energy of the batteries 1f the submarine 
wants to move for a given period [0, 7] with a speed u(t) is given by : 


EI) s Eror (Equation 4) 


Етот 15 ће total battery energy when the submarine begins to escape at time 0. 
This energy, according to assumption b), should be very close to maximum charge. ETOT 


would be normally measured in [K WHR], but for convenience we choose units for energy 


such that J= 1 in Equation 3. Therefore, we can compute Eror as U^ multiplied by the 
endurance at this speed, where U is any arbitrarily chosen speed. A realistic value for 


Eror for a modern diesel-electric submarine is: 
(16.4 mphy -(1.2 hr) = 5293 miles hrs? (Equation 5) 


C. MATHEMATICAL FORMULATION 


In Equation 1, the detection rate for a random search in an expanding disk is 
modeled as (VW)/(zy(t)). Therefore the average number of detections in the time 


interval [t, T] when the search is continuous is: 
T 


[Ea 
yy 


where a = 2RV/n. 
However our helicopter searches by dipping in an expanding circle, so the value of alpha 
will depend on the position of the dips and will be changing as a function of time. 


Therefore a better expression for aft) ıs: 


(= rateof area covered Е TR” _ К 16 


7 dl 54 20) | ER S0) 
V V 
where C is a factor near 1 that allows computing the average distance to be flown from 


dip to dip and dis the time to make a dip. If we make C=1, the average number of 


detections in the time interval [z 7] becomes: 


Г 2 il El : S Е (г); 
8 "0 бу É y(t)? e 


кә 


where K = Ди 
Ó 


A logical objective function for the submarine 1s to minimize the expected 


number of detections during the interval [ z, 7]. 


PROGRAM 1 
Minimize 


Г a(t at 
y(t) 


N 


Subject to 


|| u(tYdt € Eror and y(0) =0 


Before time т the submarine should escape with constant speed y(r)/r. Therefore the energy 


constraint for Program 1 is: 


т °?у(т)' + Ку < Етот 


This program can be expressed in terms of a Lagrangian unconstrained function as: 


F(y(t),u(t),t) = ioc + Au(t)? (Equation 7) 
y(t)” 


which is equal to 


F(y(t),u(t),t) = ЕЕ гс + Au(ty (Equation 8) 


e 
> 


Euler's equation from the calculus of variations [Ref. 8, p 399], states that a necessary 


condition for optimality is: 


F(y(t), u(t),t) — u(t) Uhren = (Equation 9) 


Euler also requires that 


с 


абе тт) UNO (Equation 10) 
Then 
L- с + Au(t) - 3Au(ty = 2 ТУ 2Au(ty; t<t<T (Едиайоп11) 
ВА N 
E ty + — 

y(t) 3y y(t) Sy 

and solving for u(t) 
1/3 
u(t) = E E - un (Equation 12) 
MA T 


If we make EUR and Roe 
2A 2A 


we have 


4 
u(t) = BY”? B с = r<ı<T (Equation 13) 
M 
V 
Á 2 1/3 
If 5 —Y and B '-v, u(t) becomes: 


Га 


ye 
u(t) = y ——————- - 1] ; rst<sT (Equation 14) 


г 
voy 4307 
Y 1s an adopted constant with units of distance and v has units of speed. 
Given y(7), Equation 14 is a first order ordinary differential equation that can be solved 
for y(t); T <t <T. 

Equation 4 states that E(T) € Eror.. For this reason the submarine’s speed might 
be any value less than or equal to u(t). Therefore, the submarine’s possible locations 
could be considered uniformly distnbuted over the circle of radius y(t) (FOC). 
According to Washburn [Ref. 5, p. 2-10], the probability of detection as a function of 
time for a random search, given y(t), can be expressed as: 

PDXty-l-e Ерот) тегт (Equation 15) 
where the random variable X is the time to detection. 


Even though Equation 14 seems not to have an analytical solution, it is possible to 


solve it numerically using the solver feature of the Excel spreadsheet. This is done by 


minimizing the value of the integral COT у(х)?)ах in Equation 15 at time 7, subject 


to the constraint imposed by Equation 4, and changing the values of y(t), v and Y. 

It is unknown whether this problem possesses the convexity properties that would 
guarantee global optimality. However, when the problem was solved repeatedly using a 
wide range of different starting values for the variables y(x), v and Y, it always converged 
to the same solutions. We suspect, but can not guarantee, this is a global optimum. 

Equation 14 might be seen as a possible speed function for the submarine and 
Equation 15 is the probability of detecting it when the search is executed inside the FOC. 
However, the helicopter doesn’t know what is the enemy is really doing, so it has to 
search covering every possible target’s behavior. The helicopter must assume that a 
possible objective function for the submarine is to maximize y(f) at every time? by 
using all of Eror. This assumption is expressed in Program 2, which has t as parameter for 


Occ Г. 


PROGRAM 2 
Maximize 
y(t) 


Subject to 
[асах < Етот 


Since the integrand is a convex function of u(x), the speed should be constant at u` (t) for 
Ü sos апа 

Ш-Н 0) T (Equation 16) 

therefore 
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пре (тома) ers (Equation 17) 

Multiplying by t 

u’ (Ot = y` (t) = Eror' 7??? 0<t<T (Equation 18) 

y'(t)is the radius of the MPD at time r. 

For any given moment f, 0 <¢ <7, u'(t) is the optimal constant speed before f to 
maximize the distance from the datum. The helicopter does not guess any particular value 


for ¢ and makes the pessimistic assumption that for every moment in the interval [0,7] the 


submarine is able to reach y’(f). 


Now let CAO (7 (г) 
dt 
Then, 
U= Бог (Equation 19) 


(Г (г) is an infeasible “Virtual Speed” function to cover every possible use of the 


submarine's battery. 


Finally, the probability of detection when the helicopter searches within the MPD 


can be expressed as 





б а(х)ах 
РР (1) =1-е б OY Lr«rteT (Equation 20) 
The following graphs show the evolution of some of the variables in Programs 1 
and 2 when considering the following realistic data: 


Constants: 
- Helicopter' s delay т = 10 minutes. 


- Sweep width W = 3 NM. 


- Submarine energy Eror- 5293 miles hours ^ . 
- Helicopter's speed V = 60 knots. 
- ö = 0.05 hours. 
- T= 2.83 hours. 
Optimized Variables: 
- Initial distance submarine-datum у(т) = 4.2385 ММ. 
- Y = 65.7236 NM. 
- v = 5.64847 knots. 
Figure 2 shows the distance y(t) from the flaming datum to the FOC that 
minimizes the probability of detection for the submarine (Program 1), and the MPD 
considered in the helicopter’s search (Program 2). Euler’s method is used to solve 


Equation 14 numerically. The distance y(t) is utilized to compute the submarine speed 


u(t). Then, the computed speed is used to obtain the distance y(t + Ar) considering that: 


y(t + At) = y(t) +u(t)At, where At =0.02 hours. (Equation 21) 
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Figure 2. Submarine distance as a function of time. 


Figure 3 shows the submarine’s velocity functions U" (rt), u(t) and u` (t). u(t) is 
the velocity function to reach the FOC shown in Figure 2 (Equation 14). U (t) is the 
submarine's virtual speed that determines the MPD assumed by the helicopter (Equation 
19). Finally, for the function u’ (t) every value is independently feasible, and each one 


implies the total consumption of the battery at time ¢ (Equation 17). 
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Figure 3. U" (t) , u (t) and u' (t) when the helicopter's delay is 
10 minutes and 7= 2.82 hrs. 


Figure 4 shows the cumulative distribution function (CDF) for the time to 
detection asa function of time when the helicopter's delay is 10 minutes and the search is 
executed: 

- Inside the MPD circle (Equation 20). 


- Inside the FOC (Equation 15). 
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Figure 4. Detection probability if helicopter’s delay is 10 minutes. 


Figure 5 shows the decrease in probability of detection as a function of the 
helicopter’s delay t, if the search is executed inside either the FOC or the MPD. This 


graph is obtained from the analytical model for the parameters listed in page 15. 
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Figure 5. Analytical Probability of Detection versus Helo”s delay. 


To obtain a more accurate result a discrete event simulation using the Java-based 
Simkit Software package is developed next, which considers the same characteristics for 
the submarine and helicopter used in the analytical model. The next chapter shows the 
methodology used to solve some of the problems related to the movement and strategies 


for the searcher and the target. 
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Ill. THE SIMULATION 


A. CONSIDERATIONS 
К Submarine Strategy 


a) Submarine Initial Position 


The initial distance from the datum y(7) is an optimized variable for each 
value of x. To determine the initial position of the submarine for each replication, the 
distance y(7) 1s multiplied by a different right triangular random varaate [0, 1] called o. 
The pair (y(r)a, 9) is then transformed to rectangular coordinates and thus creates an 


initial position uniformly distributed inside the circle of radius y(t). 


b) Velocity function 


The submarine's general velocity function is generated through 
Equation 14 


1/3 
y 
cendi 
y(t) 


D) + 2 
y(t) Sy 


u(t) — v ecc Г 


where Y andv are optimized constants for the particular helicopter's delay zand y(t) 15 
the submarine's distance from the datum at time ¢, obtained from Equations 14 and 21. A 
vector holds the velocity whose value is updated every 9 seconds in order to create a 
function which is relatively continuous. For each replication, the speed function is 
multiplied by the same c used to determine the submarine’s initial position. In doing so, 
a different speed behavior is generated for each simulation run and it is consistent with 


assumption a) for the submarine 1n Chapter II. 
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с) Course 

The submarine flees from the datum with a random course 9, restricted 
turn capability and executing 9-second legs. When a leg has finished, the submarine 
decides if it will change course. This decision is made by sampling a uniform [0,1] 
random variate and comparing it with a reference value. This reference value 1s very 
close to 1 and the maximum change in course is restricted to 3 degrees port or starboard. 


In doing so, the submarine maintains a fleeing attitude. 


2: Helicopter Strategy 


a) Position of the First Dip 


The helicopter's first dip is uniformly distributed in the circular area 
whose radius is y'(r — time to dip), where y'(t) is the MPD function obtained from 


Equation 18. 


b) Position of the Next Dip 

The helicopter dips uniformly in the area whose radius at time fis y(t’). 
Washburn and Thomas [Ref. 6,p. 420] state that, considering X as the present helicopter's 
position in the unit circle and : as the present time, the earliest time ?' to reach a point X 


in the same circle guaranteeing uniform search can be computed as: 





''s t(d - abcr e)((1- b) s T(X,X',r) (Equation 22) 

Where 

p Z А | b= с c =Cos(angle between X and Х”) 
D 2 2 2 72 212 2 

d=1+— and e =a +bd’-ab(l-c’)-2abcd 


t 


m 
ho 


U and V are the submarine and helicopter's velocity and D is the time required for the 
helicopter to make a dip with the sonar. 

The result stated in Equation 22 makes it possible to compute the position 
of the helicopter's next dip and to be certain that it will be uniformly distributed in the 
submarine's possible positions disk. Even though this procedure is applicable when 
dealing with a submarine with a constant maximum speed, a similar approach can be 


done when considering that for any time 7, the submarine’s virtual speed value U’ (t) 


may be obtained from Equation 19. 


Let 
Ages i | 


be the helicopter's position coordinates in the unit circle at time ¢ and 7 respectively, so 
x! - y! land x" -y" & 1. 


Then, the time for the next dip ¢ should satisfy 





T E) Ore шал о АОВ, 


Where 4/[x' y' (£) ^ xv (Ol -[y' y (7) - уу (91 is the distance between the helicopter’s 

present dip position at time ¢ and the future dip position at time z'. Even though Equation 

18 allows obtaining the value of y' (7), replacing it in the equation above creates a sixth 

degree equation that is not useful to solve for t'. However, we can obtain the approximate 
distance between [x, y] and [x', y'] if we assume that an approximation for y”(t') is: 

y (0) = y ()*-U (1) 


where U" (1) is the virtual speed introduced in Equation 19. 
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Replacing y' (1) in Equation 23 we have: 





V 


t+D+ 


(Equation 24) 


If both sides of the equation are multiplied by V, then 





+ + уху) + О (0е-0)-ху (ОГ + Цу OU OC-))- yw OF = Ve 
(Equation 25) 

and 

V(t* D-0) 2 -J[x0ly O+ U0 к OF +O + UV OC-D)- w’ OF 





(Equation 26) 

Squaring and reducing, then 

Va«D-ry 2pei qu qe-))- x OF tb" (04 U' O-9)- XY 
(Equation 27) 


Now let 

ttDz-a, y(t)(x-x)-xU' (tu 2b, апа y' (rn(y'-y) - y'U' (ni -c. 

Then 

V'(a-t'y [b- x'U' (QT -[co- y'U' (rp (Equation 28) 
and reducing 

Va -2aV?t«V?t? z b? x 2bx' U' (0t U" (Dy t? 9c) & 2cy'U" (0t QU" (Dy t? 


(Equation 29) 
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Now factorizing 

i? (V? - onu qy - quy )-r(gav? «2bxU* (0) - 2e! U ())- Và? - 8 -с? = 0 
а 0) 

If 

d «V^ -(xU'()y -(yU' (0, e - 2aV* - 2bx' U" (t) - 2cy U" (t) 

and f =V’a’-b’-c’ 

then, dt” -et'+f = 0 (Equation 31) 


which is a normal second degree equation. Finally 


+ е? - 44 
а СЕБ А (Equation 32) 


2d 
In practice we choose the smallest positive root of Equation 32, if greater than the 
current time ¢. The chosen time will be the solution for the earliest time that a point in the 


unit circle can be reached to guarantee a uniform dip location. 


Once :' is obtained, y' (1) can be calculated using Equation 18. Therefore, the 


coordinates of the next helicopter's dip are: 
X ny «y y y (0) 


Figure 6 shows two typical sequences of dip times when the helicopter's delay 15 


10 minutes. 
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Figure 6. Helicopter's typical dip times when delay is 10 minutes. 


3. Simulation Inputs 


a) Parameters 


In Equation 14, the parameters Y, v,and y(r) are optimized for different 
values of rt using a spreadsheet and a non-linear programming solver to solve Program 1, 
as explained in Chapter II C. The following table shows the parameters introduced in the 
simulation for each value of 7. All the columns except for the first were calculated based 
on the following constants: 
Helicopter dip radius R = 1.5 NM. 
Submarine energy Eror = 5293 miles hours" . 
Helicopter speed = 60 knots. 
Time to make a dip 6= 0.05 hours. 


At for the spreadsheet = 0.02 hours. 
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Helicopter endurance = 2.65 hours. 







PA AAA 
40] 11.48289288) 7.0162599| 101.052399 
105.7943132 
50) 13.47854091)| 7.259012 110.3977112 
114.8828 
____60] 1536430179] 7.45277030|  — 119.2756334 
124.0679351 
127.8289081 
75| 18.03435462| 7.69262138 132.0131104 


Table 1. Simulation Parameters. 


















4. Simulation Validation 
a) Average PD (PDA) 
The average probability of detection (PD) for a particular value of the 
helicopter’s delay rıs obtained if it is considered that for a helicopter dip at time г, Ше 


probability of detection might be computed as: 


4 


- 


уа) 





PD(i" dip) = (Equation 33) 


Let N be the number of dips made by the helicopter, including the detecting dip (1f any) 
but none beyond it. Then the probability of detection PD for any simulation replication 


1S: 
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N 
BI 29) PD ap») (Equation 34) 


ігі! 
The expected value is estimated by replication. Letting V, be the number of dips in the 
k" replication and ¢,, be the time of the 7" dip in the k™ replication, the estimate for Q 


replications 15: 


о М, R 





Average PD (7,0) = аи (Equation 35) 


B. BASIC RESULTS 


10000 replications were run for 14 different values of т. Тһе simulation 
simultaneously computed Monte Carlo estimates for the value of PD, . Figure 7 (next 
page) shows the evolution of the probability of detection and its 95 % confidence 
intervals versus tas estimated by the simulation compared with PDa, when the 
helicopter’s search 1s executed inside the MPD. It can be seen that the values of the 
estimates PDa remain inside the confidence intervals for all values of t, confirming the 
accuracy of the simulation results. In the next chapter the simulation is used to predict the 
probability of detection when the behavior of the submarine and the helicopter are varied, 
for helicopter’s delay of 10, 20 and 30 minutes. 

Tables with detailed simulation results used in this chapter are included in 


Appendix A. 
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Figure 7. Probability of Detection obtained from simulation and estimate 
PDA versus helicopter's delay, when search is inside MPD. 
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IV. RESULTS AND ANALYSIS 


А. MODEL FITTING 


Figure 8 shows the probability of detection as a function of the helicopter's delay 
t according to the results obtained from the simulation, when the search 1s executed 
inside the MPD. A fitted function is also shown and the corresponding regression 


analysis is detailed in Appendix B. That curve is valid for the parameters listed in chapter 
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Figure 8. Simulation results and fitted function. 
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B. TACTICAL IMPLICATIONS 
Different variations to the problems were explored. The following results were 
obtained from 10000 simulation replications. A large-sample 100(1-a)% confidence 


interval for a population proportion p, if n is the sample size is: 


^ ^ 


MT Dum) 


m2 
n 


^ 


p 


^ ^ 


When we don't know the value of p the most conservative approach is to use p =.5, to 


maximize the value of the expression under the root [Ref. 9, p.289]. Therefore a 9596 


confidence interval is: 
Е: e — = +0.01 
10000 


р Submarine on the Edge of the FOC 

The submarine might decide to minimize the analytical probability of being 
detected by remaining on the edge of the FOC. Figure 9 shows the probability of 
detection for helicopter delays of 10, 20 and 30 minutes when the submarine escapes by 
sailing on the edge of the FOC compared with when its position 1s uniformly distributed 
in the same area. 

A small but statistically significant decrease in the probability of detection is 


observed for helicopter's delay of 10 minutes (p value = 0.000975). 
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Figure 9. Submarine escaping on the edge of the FOC compared 
with when its position is uniformly distributed in the same area. 


This effect seems to be attractive for the submarine. However, according to our model, 
this tactic has the disadvantage that at the end of the helicopter’s search period the 
submarine's battery 1s completely exhausted and its options are reduced to being 
immobile, at the bottom or snorkeling. No statistically significant differences are 
observed for delays of 20 and 30 minutes (p values = 0.51 and 0.15 respectively). The 
tactical significance is small in all cases, which should not be surprising because the 


helicopter makes no attempt to exploit the predictability of the range to the submarine. 
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2. Submarine Remains in the Vicinity of the Datum 

The submarine might attempt to surprise the enemy force by remaining in the 
vicinity of the datum. This situation may be also presented if the submarine executes an 
involuntary attack and is forced to remain close to the datum because of its low battery 
level. Figure 10 shows the effectiveness of this tactic for 3 different values of the 
helicopter's delay compared with when the submarine behavior determines a uniformly 
distributed position in the FOC. If the submarine remains in the datum vicinity, no 
statistically significant differences are observed for delays of 10, 20 or 30 minutes (p 


values = 0.31, 0.43 and 0.95 respectively). 
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Figure 10. Probability of detection when the submarine remains 
in the datum vicinity, compared with when its position 1s uniform. 





These results indicate that the submarine does not get any tactical benefit 
remaining in the vicinity of the datum after the attack. Again, this should not be 


surprising given the helicopter's tactics. 
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3. Helicopter Fails First Detection Success 

When the helicopter is dipping, a poorly trained crew may fail to detect the target 
even when the target is within range. When dealing with a diesel-electric submarine this 
is especially possible. They almost make no noise and most of them are small. Although 
it is logical to suppose that this lack of skills will negatively affect the probability of 
detection, our intention is to measure its real impact in the mission accomplishment. In 
Figure 11, a dramatic loss of probability of detection in the event that the helicopter fails 


to detect the first time the submarine is detectable, for all helicopter’s delays. 
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Figure 11. Probability of Detection when helicopter fails to detect 
the submarine when an opportunity 1s presented. 


In this case it 1s very unlikely that the helicopter will ever have a second chance of 


detection. 
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4. Helicopter Searches Inside the FOC 

The helicopter’s search may be executed inside the circle determined by the FOC 
instead of the MPD. In that case, the helicopter assumes correctly that the submarine will 
remain inside the FOC the whole time and the search is conducted according to this 
assumption. When this assumption Is correct, a small increase in probability of detection (p 


value z 0) 1s observed for all three delays (Figure 12). 
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Figure 12. Probability of detection when the helicopter's search is inside 
the MPD or the FOC. 


However, assuming that the submarine is inside the FOC will be incorrect whenever the 


submarine exceeds the limits of that circle. To do that the submarine could follow any of 


many speed behaviors determined by u'"(t) given by Equation 17. For these cases 
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neither the analytical model nor the simulation allow us to know the effects in terms of 


the probability of detection of searching based on a wrong assumption. 
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V. DISCUSSIONS OF THE MODELS AND FUTURE WORK 


A. ANALYTICAL MODEL 


The analytical model explored in this thesis has some evident weaknesses: 


le Capacity of the Submarine Battery 
In Chapter II it was stated that the constraint imposed by the battery's’ energy if 
the submarine wants to move for a given period [0, 7] with a speed u(t) is: 


ВИИ TOT (Equation 4) 


Eror is computed as U^ multiplied by the endurance at this speed, where U is any 
arbitrarily chosen speed. For an ideal battery this should be true for every chosen U. 
However, batteries do not always behave in an ideal manner and therefore this method 
of obtaining the value of Eror 1s not exact. It 1s suggested that future studies may focus 


on obtaining a better model for the battery. 


2. The Edge Effect 

The analytical model explained in Chapter II fails to capture the effect produced 
when the helicopter dips in a position where the distance to the datum S is greater than 
the submarine FOC minus the sonar range R. When this is the case, the detection rate 
decreases because the area of the sonar detection circle A(S) inside the submarine 
possible position disk also decreases (Figure 13). This effect is a source of inaccuracy for 
the model and for the probability PD explained as a simulation verification method in 
Chapter III. However, it does not seem to be very significant. Future studies may try to 


measure its effect in terms of probability of detection and confirm this belief. 
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Figure 13. The Edge Effect is produced when the dip radius S is outside y(t)— R. 


3. The First Dip 

The helicopter's first dip has a significant value in terms of probability of 
detection. The models studied in Chapter II gradually integrate the probability of 
detection during the interval [z, 7], where r is the helicopter delay and T is the ending 
time of the search period. Thus, they do not capture properly that an important part of the 
probability is added almost immediately during the first dip. This is a source of difference 
between the analytic models and the simulation. Table 2 shows the probability of 
detection observed from 10000 replications of the simulation for different helicopter 


delays and the magnitude of such a probability coming from the first dip detection when 
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the search is conducted inside the MPD circle. It is evident that a significant percentage 


of the total probability of detection is because of the first dip. 


[Helo Delay | Alldips | Firstaip | 
40 f 00355 | 0.0124 | 
a | oe ШЕТ ЖЕ 
Table 2. Helicopter's First Dip. 





B. FURTHER STUDIES 

This study represents a more realistic method of posing the Flaming Datum 
problem when dealing with diesel-electric submarines that must respect an energy 
constraint. This constraint was never considered in any other study. For this reason the 
developed models may be considered a valid starting point for future works. It is 
suggested that they should be focused in the following aspects: 

Incorporating the first dip effect as an important source of probability of detection 
for the analytical model. 

Improve the modeling of the submarine battery. By doing so, more accurate limits 
for the submarine capabilities will be established and the credibility of the predictions of 


the two models will be increased. 
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APPENDIX A. TABLE OF SIMULATION RESULTS 


The table shows the probability of detection as a function of the helicopter's 
delay. The second column corresponds to the values obtained from the analytical model 
when the helicopter's dips inside the MPD area. The third, fourth and fifth columns show 
the values generated by the simulation and its 95 % confidence intervals. The last column 


is PDa, which was used to verify the credibility of the data generated by the simulation. 
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APPENDIX B. REGRESSION FOR THE PROBABILITY OF DETECTION 


The following output shows the regression model for the probability of detection 
as a function of the helicopter’s delay t according to the results obtained from the 


simulation, when the search is executed inside the MPD. 


*** Linear Model *** 


Call: Im(formula - probability - delay^3 « delay^2 + delay, dat 
regression.frame, na.action 
- na.exclude) 


Residuals: 
Min TO Median 30 Max 
-0.00502 -0.001252 0.0002859 0.001655 0.005129 
Coefficients: 
Value Std. Error t value Pr(>|t|) 
(Intercept) 0.2044 0.0074 27.5633 0.0000 
I(delay^3) 0.0000 0.0000 -6.4155 0.0001 
I(delay^2) 0.0002 0.0000 8.5822 0.0000 
delay -0.0090 0 0007 = 23. 2313 0.0000 


Residual standard error: 0.003216 on 10 degrees of freedom 
Multiple R-Squared: 0.9935 

F-statistic: 506.2 on 3 and 10 degrees of freedom, the p-value 
3.234е-011 


Analysis of Variance Table 
Response: probability 


Terms added sequentially (first to last) 
Df Sum of Sq Mean Sq F Value Pr(F) 
I(delay^3) 1 0.007579049 0.007579049 732.5684 1.096000e-010 
I(delay^2) 1 0:006321668 0.0063216693 611.0536 2.6820006-010 
delay. 1 0 001811225 0:001811225 17520676 1. 160015е-007 
Residuals 10 0.000103459 0.000010346 
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